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Abstract 

We explore imprints of the solar wind interaction with the local interstellar medium 
on spectral properties of the backscattered solar Ly-a radiation. We employ newly 
developed effective model for the interstellar H atom velocity distribution function 
in the heliosphere (Katushkina & Izmodenov 2010, 2011). This model takes into 
account both global effects of H atom perturbations at the heliospheric boundaries 
and local (i.e. within 10-20 AU from the Sun) effects of the solar ionization, charge 
exchange, solar gravitation and radiation pressure. 

Backscattered solar Ly-a profiles and their zero, first and second moments were 
computed at 1 AU for the anti-solar directions of line-of-sight. Then the moments 
are compared with the moments calculated by using more simplified one-component 
and two-component hot models. The comparison shows that the Ly-a intensities 
calculated in the frame of new model are somewhat smaller than those calculated 
by means of the simplified models. Small differences in the first moment (i.e. line- 
shift) are also observed in the downwind direction. 

The largest difference between new model and hot models is seen in the second 
moment of the backscattered Ly-a profile. This moment is also called as line-of- 
sight "temperature" or line-width of the backscattered spectra. The spectral width 
depends on the width of the H atom velocity distribution function within few AU 
from the Sun. We show that the line- widths calculated on the basis of the our model 
is qualitatively different from those of simplified models. 

It is also shown that minimum of the line-width at 60° from upwind that was 
observed by SOHO/SWAN (Costa et al., 1999; Quemerais et al., 2006) is absent 
in the frame of our model. Such a minimum is expected from simple consideration 
of two populations (primary and secondary) of interstellar H atoms and although 
new model includes two populations of H atoms we do not see this minimum in 
our results. Instead of this minimum, our model predicts a small local maximum of 
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the line-width at ~150° from upwind. We explain these effects by non-Maxwellian 
behavior of the velocity distribution function of H atoms at entrance to the helio- 
sphere. More specifically, difference in the kinetic temperatures T x and T z plays a 
key role. Physical phenomena which can help to get the observed minimum in the 
frame of our model are discussed. 

Key words: heliosphere, interstellar hydrogen, hot model, solar Ly-a radiation 
PACS: 96.50.Xy, 96.50.Zc, 98.70.Vc 



1 Introduction 



Our solar system is moving though the Circum-Heliospheric Interstellar 
Medium (CHISM). The solar wind interacts with interstellar plasma and com- 
plex structure of the interaction region is formed due to relative Sun/CHISM 
motion (see FigJT] A). CHISM is a weakly (20-30 %) ionized plasma and 
mainly consists of hydrogen. The mean free path of the interstellar hydro- 
gen atoms (H atoms) is larger or comparable with the characteristic distance 
of the SW/CHISM interaction region. Therefore, H atoms penetrate deeply 
to the heliosphere, where properties of H atom component can be measured 
directly on indirectly. Due to their large mean free path in the heliosphere a 
kinetic approach should be used to describe H atom component theoretically. 

The first self-consistent kinet ic-gasdynamical model of t he SW/CHISM in- 
teraction has been developed by Baranov Sz Malamal f 1993h (Baranov-Malama 
model hereafter). It was shown that the SW/CHISM interaction region (that 
is often called in literature as the heliospheric interface) consists of four regions 
with different plasma properties. These regions are separated by three discon- 
tinuity surfaces, namely, the termination shock and bow shock are boundaries 
of the supersonic solar wind and interstellar medium flows, accordingly, and 
the heliopause is a tangential surface separated charged components of the 
two plasma flows. 



One of the results ob tained by iBaranov &: Malamal ( 119931 ) and previously 
by IBaranov et al.l (1199 lh is the hydrogen wall that was discovered experimen- 
tallv bv lLinskv fc Wood! (1l996h (see also Wood et al, 2007a, 2007b for recent 
developments). The hydrogen wall is an increase of H atom density in the vicin- 
ity of the heliopause due to production of newly created H atoms by charge 
exchange. These secondary atoms have properties of plasma in the vicinity 
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Figure 1. A: Qualitative picture of the heliospheric interface; B: coordinate system 
used in the models: z is the axis of symmetry directed toward the upwind direction, 
p is an axis in cylindrical coordinate system that is perpendicular to the upwind 
direction; curve 1 is the outer boundary in our model; curve 2 is a schematic repre- 
sentation of the heliospheric termination shock. 

of the heliopause. The interstellar plasma in the vicinity of the heliopause is 
heated and decelerated, and, therefore, the mean velocity of the secondary H 
atom population is smaller and the temperature is larger as compared with 
the velocity and temperature of the original (or primary) H atom population 
that consists of interstellar atoms, which entering into the heliopshere without 
being charge exchanged in the heliospheric interface. A mixture of the primary 
and secondary interstellar H atom populations penetrate to the heliosphere. 
Inside the heliosphere in region between the heliopause and the termination 
shock (this region is often called as the inner heliosheath) new atoms are 
also created during charge-exchange. These atoms are called as the energetic 
neutral atoms (ENA s) and now they are m easured by Interstellar Boundary 



Explorer spacecraft (IMcComas et al.l . 120091 ) . ENAs have a small number den- 



sity in the heliosphere as compared with interstellar atoms, that is why ENAs 
are not important for backscattered Ly-a radiation. 

It was recognized since 1970s that the solar Ly-a photons are effectively 
backscattered by the interstellar H atoms. Spectral properties of the backscat- 
tered Ly-a radiation depend on the velocity and spatial distribution of the 
interstellar hydrogen inside the heliosphere. Since the distribution of H atoms 
is disturbed in the heliospheric interface, then the backscattered Ly-a may, 
in principle, serve as t he source of i nform ation on the heliospheric interface 



structure. For example. ICosta et al.l (119991 ) analyzed SOHO/SWAN data and 
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showed that the temperature of the interstellar hydrogen at large distances 
from the Sun (at 50 AU where local effects of solar gravitation, radiation 
and ionization are negligible) is significantly higher then temperature of the 
CHISM. The bulk velocity of the hydrogen flow obtained by this analysis is 
about 21 km/s that is ~5 km/s smaller than the relative SW/CHISM ve- 
locity. These differences in the properties of the hydrogen gas as compared 
with the properties of the CHISM are explained by existence of the secondary 
interstellar H atom population in the heliopshere. 



The effect of the heliospheric interface is also seen in the measured ( ILallement et al. , 



20051 ) deflection of the interstellar H atom flow direction as compared with the 
direction of the interstellar helium. In contrast to the hydrogen the charge- 
exchange cross-section of helium {He + H + <-»• He + + H) is small, that is why 
it is usually assumed that the interstellar helium stay practically undisturbed 
during its motion from CHISM to the heliosphere. It means that the velocity 
vector of the helium flow in the vicinity of the Sun refers to the value and 
to direction of the interstellar gas flow. In turn, hydrogen atom flow is af- 
fected by charge-exchange with protons in the SW/CHISM interaction region. 
The deflection of hydrogen atoms inside the heliosphere from the direction of 
SW/CHISM relative motion is most probably due to distortion of the global 
heliospheric interface structure caused by the interstellar magnetic field that 
leads to asymmetry of interstellar plasma flow in the heliospheric interface and 
this asymmetry transfers to H atoms through the charge-exchange. This fact 
opened a possibilit y to estimate magnitud e and direction of the interstellar 
magnetic field (e.g. Ilzmodenov et al.l . 120051 ) . 



Another manifestation of the heliospheric interface was found during ana- 
lyzes of the second moment of the backscattered Ly-a radiation spectra. This 
moment corresponds to the width of the spectra and depends on the width 
of the H atom velocity distribution. Being expressed in the units of temper- 
ature the line-widt h is often called as th e line-o f-sight (or 1-o-s) temperature. 
Costa et al. ( 1999 ) and Quemerais et al.l ( 2006 ) derived the 1-o-s temperature 
of the backscattered Ly-o; from SOHO/SWAN data and presented the temper- 
ature as a function of the angle 9 between the line-of-sight direction and the 
direction toward the CHISM flow (i.e. the upwind direction). It was shown that 
the temperature in the directions of 9 = 50-60° is smaller by 1500- 2000 K than 



the te mper ature toward upw i nd. Th is minimum was explained by ICosta et al 



( 119991 ) and I Quemerais et al.l ( 120061 ) qualitatively by presence of the primary 
and secondary populations of the interstellar hydrogen atoms. 



The goal of this study is to theoretically explore the imprints of the he- 
liospheric interface on spectral properties of the backscattered Ly-a radiation 
and separate these effects from the local effects (i.e. the effects that act within 
10-20 AU of the Sun) of the solar ionization, charge exchange, solar gravitation 
and radiation pressure. We have computed line profiles and their moments of 
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the backscattered Ly-a spectra based on newly developed advanced kinetic 
model of the hydrogen distribution in the heliosphere (Katushkina & Izmode- 
nov, 2010, 2011). This model combines comparative simplicity of the classical 
hot model with ability to account for both the perturbations of the hydro- 
gen distribution function in the SW/CHISM interaction region (later it will 
be mentioned as effects of the heliospheric interface) and local effects which 
are significant only in the vicinity of the Sun. In order to separate the helio- 
spheric boundary effects from the local effects we compare new model results 
with the results of the simplified hot models, which describe local effects in 
similar way as the advanced model, but do not take into account the effects 
of the heliospheric boundaries. 



2 Modeling of the hydrogen distribution inside the heliosphere 



To model spatial and velocity distribution of the interstellar hydrogen in 
the heliosphere we solve a kinetic equation (see, e.g. Katushkina & Izmodenov, 
2010) for the distribution function /(r, w, t) of the hydrogen atoms inside the 
sphere with radius of 90 AU around the Sun (see Fig. [1] B). The processes 
of the photoionization, charge-exchange, solar gravitation and solar radiation 
pressure are taken into account. It is assumed in the model that the rates of 
charge exchange (3 ex and photoionization (3 p h decrease with distance from the 
Sun as 1/r 2 . We use following values of ionization rates at the Earth orbit: 
P P h,E = 116 x 10~ 7 s" 1 , (5 £X) e = 4.8 x 10~ 7 s _1 . To calculate the charge- 
exchange rate at the Earth orbit we use averaged values of the solar wind 
number density (6 cm -3 ) and velocity ( 440 km/s) known from mea surements, 



and charge-exchange cross-section from lLindsay fc Stebbingi (120051 ). The solar 



gravitational force and the force of the radiation pressure act in opposite 
directions and both are proportional to 1/r 2 , where r is a distance to the Sun. 
Dimensionless parameter /i defines a balance between the solar gravitation 
and radiation pressure. In this work we assume that [i is equal to 1.258. 

In this paper we consider three models which differ each from other by 
boundary conditions at 90 AU. Model 1 is the one-component hot model 
that implies the Maxwellian velocity distribution function at the outer bound- 
ary. For the specific calculations performed in the frame of this model we use 
parameters of the Maxwellian distribution (i.e. number density, bulk velocity 
and temperature) calculated at 90 AU in the upwind direction in the frame 
of the Baranov-Malama model with parameters of chosen to be the same as 
in Katushkina fc Izmodenovl (120101 ) . These values at 90 AU are: n# = 0.54 
■n H ,LiSM, Vh, z = -0.79-V h ,lism, Th = 15 000 K, where n H ,LiSM=0.18 cm -3 , 
V# i ljsm=26.4 km/s. These parameters correspond to the mixture of the pri- 
mary and secondary interstellar atoms. Note, the effects of the heliospheric 
interface are taken in this model at zero level by taking parameters in the 
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Maxwellian distribution different from the CHISM parameters. 



Model 2 is so-called two-component hot model (e.g. iBzowski et all 120081 ) . 
In this case model 1 is employed separately for the primary and secondary pop- 
ulations of interstellar atoms. The velocity distribution functions are assumed 
to be Maxwellian for each of the populations and parameters of the distribu- 
tion functions were calculated in the frame of Baranov-Malama model. These 
parameters for the primary interstellar atoms are: n Hprim = 0.22 -n HjLISM , 
V z ,prim = -1-06-Vh,lism, ^H^rim = 6 840 K, and for the secondary interstel- 
lar atoms: n H ^ second = 0.32 -tih^lism-, "V \, second = -®-6&-Vh,lismi T# jSeconrf = 18 
126 K. 



As it was shown in iKatushkina fc Izmodenovl (120101 ). model 2 leads to sig- 
nificant discrepancies in the distribution of H atoms inside the heliosphere as 
compared with the self-consistent Baranov-Malama model. This fact is con- 
nected with non-Maxwellian behavior of the velocity distrib ution function of 



hydrogen atoms after th eir crossing the heliospheric interface (llzmodenov et al 
200ll : llzmodenovl . l200lh . 



In model 3 we use the boundary conditions allowing to get the distribution 
of H atom inside the heliosphere in very good agreement with the distribution 
obtained in the frame of Baranov-Malama model. For the primary interstellar 
atoms the distribution function at 90 AU is chosen as 3D normal distribu- 
tion. In this case all zero, first and seco nd moments of the velocity dist ribu- 



Katushkina fc Izmodenov . 2010l ) and 



tion function are taken into account (see 
calculated in the frame of Baranov-Malama model. However, for the popula- 
tion of the secondary interstellar atoms such approach does not give the full 
agreement with the Baranov-Malama model because the distribution func- 
tion of this population is asymmetric with respect to its maximum and it has 
nonzero third moments which are neglected in th e 3D normal distribution. 
That is why for the secondary interstellar atoms IKatushkina &: Izmodenov 
(120 111 ) have used the velocity distribution function of at 90 AU calculated by 
Monte-Carlo in the frame of the Baranov-M alama model. It was shown that 
15-20 millions of unsplitted ( iMalamal . 119911 ) trajectories in the Monte-Carlo 
code is en ough in order to get acceptable ac curacy of calculations. In this paper 
we follow IKatushkina fc Izmodenov! (120111 ) approach for the primary interstel- 
lar atoms and use the numerically calculated velocity distribution function as 
the boundary conditions at 90 AU for the secondary interstellar atoms. 



It is important to note here that for the stationary and axisymmetric case 
model 3 provides the same distributions within 90 AU from the Sun as the 
Baranov-Malama model. However, computationally our model is more effec- 
tive and allows to calculate the velocity distribution function everywhere inside 
the heliosphere with required precision. Moreover, our model allows to study 
local (i.e. within the heliosphere) time-dependent and 3D phenomena very 
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effectively. 



3 Modeling of the backscattered solar Ly-o profiles 



The backscattered Ly-o profiles I(r,z/,f2) were computed for anti-solar ra- 
dial directions (f2) for an observer located at 1 AU. Here r is a position of 
observer, v is a frequency of the backscatt ered radiation, ft i s the line-of-sight. 



We use "self-absorption" approximation (jQuemeraisl . 120001 ) . In this approxi- 
mation only once scattered photons are considered (i.e. multiple scattering of 
photons are neglected), and absorption of the photons between the Sun and 
the scattering point is neglected. Despite this simplified "self-absorption" ap- 
proach gives similar resul t s as c ompared with a full radiative transfer model. 



Quemerais &: Izmodenovl ( 120021 ) have shown that for the line-width of the 
backscattered profile at 1 AU the difference between full radiative transfer 
model and the self-absorption model is less than 15 % for the upwind direc- 
tion and becomes smaller as the line-of-sight moves from upwind to downwind. 
In the downwind direction the two approaches give practically the same line- 
widths. The simplified self-absorption approach is sufficient for the purposes 
of this paper. As it will be seen the effects of the heliospheric interface are 
essentially larger than the difference between the full radiative transfer and 
self-absorption models. 

The radiative transfer equation for I(r,i/,f2) can be written as follows: 

n ■ V/(r, v, n) = -<7„(r, u)n H {r)I{r, v, O) + n H {r)j{r, v, -O). (1) 

Here n#(r) is the number density of hydrogen atoms, cr v (r, v) is the differential 
scattering cross-section that is proportional to the projection of the hydrogen 
distribution function on the line-of-sight, j(r,i/,-Q) is the atomic emission coef- 
ficient which measures the number of photons emitted by a hydrogen atom per 
second per frequency unit and per solid angle. Note, that scattered photons 
travel in the direction opposite to the line-of-sight direction i.e. in — Q. The 
first term in the right hand of the equation ([T]) is a loss term due to absorption 
of the scattered photons and the second term is a local source of the emission 
due to scattering process. 

The equation ([!]) has a classical solution in case of the self-absorption ap- 
proximation: 

A TT /"oo 

I(r,v,n) = — y o n H (r + sn)j(r + sn,u,-n)e- T ^ +sn ^ds, (2) 

where, s is a coordinate along the line-of-sight, T u (r', r) is the optical thickness 
for scattered photons with frequency v going from the scattered point r' = 
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r + sf2 to observer at point r. Ly-a profile is measured here in Rayleigh. In this 
solution only once scattered photons are taken into account, and absorption 
of primary solar photons between the Sun and scattering point is neglected. 
Atomic emission coefficient can be represented as: 

j(r', u, -O) = F s (r', u p ) a v (v' , u) . (3) 

Here <p(oj) is a scattering phase function that expresses the relat ion between di- 



rectio n of propagation of the photon before and after scattering (IBrandt fc Chamberlain! 



19591 ). Fs(r', u p ) is a flux of primary (solar) Ly-a photons with the fre quency 



of Up at point r'. We use the solar Ly-a spectra (ILemaire et all Il998l ) to cal- 
culate solar Ly-a flux at the Earth orbit with defined frequency. In that case 
when line-of-sight is radial, there is a simple relation between frequency of 
the primary solar photon u p and frequency of the backscattered photon v: 
Up = 2 ■ u Q — u, u is the frequency at line center. Thus, if we know the hydro- 
gen velocity distribution function in the whole heliosphere, we can calculate 
the profiles of the backscattered Ly-a radiation. 

Since we consider the two different populations of interstellar H atoms in- 
side the heliosphere it is possible to calculate profiles of the radiation scattered 
by each populations separately. To do this we consider photons that were scat- 
tered by the primary and secondary interstellar atoms independently. Optical 
thickness is calculated for the mixture of the primary and secondary atoms, 
since a photon scattered, for example, on the primary interstellar atom may 
then be absorbed by both the primary or secondary atoms. 

We calculate the following moments of the backscattered radiation profiles: 



// os (r, fl) = / J(r, u, fl) du - intensity in Rayleigh; 
Jo 

Vi os (r, fl) = — — j ' ' — line-shift expressed in [km/s]; 

TUr, O) = ^ ^Hu)-VUrM)?I^uM)du ^ ^ 

Here, u(u) = c(u/u — 1), ran is the mass of a hydrogen atom, kb is the 
Boltzmann constant. The line-shift of the spectra is often called as the line-of 
sight velocity, and the line-width is also called as the line-of-sight (or apparent) 
temperature. These integral characteristics of the backscattered Ly-a profile 
reflect the properties of the velocity distribution function of H atoms inside 
the heliosphere, but they do not coincide and should not be mixed up with 
the bulk velocity and temperature of hydrogen atoms far away from the Sun. 
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4 Results of the modeling 

We computed the backscattered Ly-a profiles and their moments by making 
use of three described above models of the hydrogen distribution inside the 
heliosphere. All computations were performed for the anti-solar directions. For 
the 2D axisymmetric problem considered here each line of sight is characterized 
by one angle 9 that is counted from z-axis (see FigJT]B). 

Fig. H] shows the intensities (left column), line-shifts (center-column) and 
line-widths (right column) of the backscattered Ly-a radiation at 1 AU as 
functions of the line-of-sight angle 9. For models 2 and 3 the profiles of the 
photons scattered by the primary and secondary interstellar H atoms were 
computed separately (plots B, C in Fig. [2]). The total backscattered profile 
was calculated as well (plot A). For model 1 that has only one component of 
hydrogen we compute only characteristics of the total Ly-a radiation. 

For the total radiation (plots A) models 1 and 2 lead to systematic increase 
of the intensities as compared with the intensities calculated for model 3. 
Comparison of models 2 and 3 for the primary and secondary populations 
(plots B and C) shows that the increase in intensities is due to the secondary 
H atom component (compare left-columns in the plots A and C). Models 2 
and 3 are nicely agreed for the primary H atom component . 

Now, we will explore the main reason of the different results produced by 
models 2 and 3. The models are only different in the boundary condition. We 
can identify two differences of the boundary conditions. The first is the de- 
pendence of the boundary conditions on the angle 9 in model 3 that does not 
exist in model 2, since model 2 assumes constant parameters at the boundary. 
The second difference is the non-Maxwellian velocity distribution function at 
90 AU in model 3, while in model 2 the velocity distribution is Maxwellian. 
To explore which of the two identified reasons is mainly responsible for the 
differences in the spectral properties of Ly-a we perform additional model cal- 
culation (model 2a). In model 2a we assume that velocity distribution func- 
tions at the outer boundary of 90 AU are Maxwellain for both the primary 
and the secondary populations of H atoms. However, the parameters of the 
Maxwellians (i.e. number density, bulk velocity and temperature) are functions 
of the angle 9 that were calculated in the frame of Baranov-Malama model. 
Therefore, model 2a allows to separate the effects of non-uniform flow of H 
atoms at 90 AU after they passed the heliospheric interface from the kinetic 
effects connected with non-Maxwellian features of the velocity distribution 
function. 

Results of model 2a are also shown in Fig. [2] A. It is seen that model 2a 
and model 3 produce very close results in the backscattered Ly-a intensities, 
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Figure 2. Intensities (left column), line-shifts (center-column) and line- widths (right 
column) of the backscattered Ly-a radiation at 1 AU as functions of the line-of-sight 
angle 9. Plots A (top row) are for the total radiation scattered by both primary and 
secondary interstellar atoms. Plots B (middle row) are for photons that were scat- 
tered by the primary interstellar atoms. Plots C (bottom row) correspond to the 
radiation scattered by the secondary interstellar atoms. Different curves correspond 
to three models of the hydrogen distributions in the heliosphere: (1) is the one-com- 
ponent hot model (model 1); (2) is the two-component hot model (model 2) ; (3) 
is our model (model 3) that takes into account effects of the heliospheric interface; 
additional curves marked as 2a in plots A,B,C correspond to model of H atoms 
described in this section, it is model 2a that is equal to model 2 plus ^-dependence 
of hydrogen parameters at the outer boundary. 



although there is a difference of 20 rayleigh for the upwind direction. It means 
that the main difference of models 2 and 3 is due to non-uniform flow of H 
atoms at 90 AU in model 3. However, still there is a smaller difference due 
to non-Maxwellian behavior of the velocity distribution function at 90 AU 
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(which is taken into account only in model 3). 



At first view, it is not evident why the angle dependence of the H atom 
parameters at 90 AU would strongly influence the backscattered Ly-a emission 
measured at 1 AU. It is especially so, because the main emissivity region for 
the backscattered Ly-a radiation at 1 AU is located approximately at 2 AU for 
the upwind line-of-sight and at 7 AU for downwind. From the simple (naive) 
consideration one could expect that the most of H atoms in the regions close 
to the Sun would arrive from upwind. However, Lallement fc Bertaux ( 1990h 



have shown (in the case of the hot model) that the most of the atoms arrive 
into the vicinity of the Sun not from upwind. The same is true for our model. 
To illustrate this we calculate function n(r$, 9b) for two points (i=l,2) inside 
the maximum emissivity region. Point 1 is located in the upwind direction at 
2 AU from the Sun. Point 2 is located in the downwind direction at 7 AU 
from the Sun. The function n(r^,0&) represents the contribution to the total 
number density at a given point (point 1 for curve 1 in Fig. H] A and point 2 
for curve 2) from those particles that cross the outer boundary of 90 AU at 
9 = 9b- This function n(i"j, 9b) is defined as follows (for certainty we present 
all formulas only for point 1): 

n(ri,0 6 ) = / /(ri,wi)dwi (4) 

Here / is the velocity distribution function of H atom. The integration is 
performed over those wi that correspond to the trajectories crossing the outer 
sphere of 90 AU at 9 = 9b- Here for simplicity we consider a balance between 
the solar gravitation and radiation pressure (// = 1). This means that all 
atom's trajectories are straight lines. In this case subspace Qi is a cone with 
apex angle equal to 9\ as it is seen from Fig. [3J Let us introduce a spherical 
coordinate system in the velocity space. It means that we describe velocity 
vector w by its module w and two spherical angles 9 and <f, i.e. orthogonal 
coordinates of vector w can be represented as follows: 

w x = w ■ sin{9) ■ cos(<f) 

w y = w ■ sin(9) ■ sin(<f) (5) 
w z = —w ■ cos (9) 



In spherical coordinates dwi = w\ sin(9) dw\ d9 dip. For velocities from the 
subspace 9 — 9\ — const for chosen value of 9b and integration over 9 is 
not needed. Hence, for our case integral (j3j) can be rewritten in the following 
form: 

/*+oo r2ir _ 

n(r 1 ,9 b )= / f(r 1 ,wi,9 = 9 1 ,f)w\sin{9 l )dw l df. (6) 

Jo Jo 
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Figure 3. Schematic picture of penetration of the hydrogen atoms from the outer 
boundary to the vicinity of the Sun. Point 1 (ri) is located in the upwind direction 
at 2 AU from the Sun. In the case of straight atom's trajectories we consider only 
those atoms that arrive to point 1 from the outer boundary at 9 b . Angle 9\ can be 
determined if we know r\ and 9 b . 

Now, in the case of [i = 1: 

f(rx, wi, 9i, 0) = f b (9 b , wt, 9 1 ,(p) ■ exp(-A(r 1 , 9 b , w{)), 

where f b is a corresponding velocity distribution function at 90 AU, A is the 
loss of hydrogen atoms along its trajectory from the outer boundary to point 
1 due to ionization processes and 9i = 6\{6 b ). f b does not depend on angle (p 
due to axial symmetry of the boundary conditions. Therefore expression ([6l) 
can be represented as: 

/*+oo r2ir _ 

n(r 1 ,9 b )= / f b (9 b ,Wi,9i(9 b )) exp(-A(r 1 ,9 b ,w 1 ))wlsin(9 1 )dw 1 d<fi = 

Jo Jo 

r+oo 

= 2ix sin{6i) ■ / f b (9 b ,wi,9i(9 b )) exp(-^l(ri, 9 b , wi)) w 1 dwi = 
Jo 

= 2nsin(9 1 ) ■ ^(ri,^). 

Function f b (9 b ,Wi) decreases with 9 b for each given value of Wi, because 
at 90 AU V Zt H ^ V Pt H and distribution function has a maximum in the up- 
wind direction. Loss-function A(r 1; 9 b , w^) increases with 9 b for point 1 because 
length of the trajectory has a minimal value in upwind and exp(— A) decreases 
with 9 b . Hence, function g(ri, 9 b ) decreases with 9 b . And sin(9i) increases with 
9 b for 9 b G [0, 7r/2]. Therefore n(ri, 9 b ) that is a product of sin(9i) and g(ri, 9 b ) 
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Figure 4. Contributions to the total number density (A) and radial velocity (B) 
of atoms arrived from various directions. Curve 1 in plot A corresponds to point 1 
(r = 2 AU, 6=0°), and curve 2 corresponds point 2 (r = 7 AU, 0=180°). Contribu- 
tion to the radial velocity was calculated only for point 2. Different curves in plot B 
correspond to different models of hydrogen distribution: solid line is model 3, dashed 
line is model 2 and dashed-dot line is model 1. All calculations are performed for 
the mixture of primary and secondary interstellar atoms. 

should have a maximum for some 9b between 0° and 90°. 

Fig.HJA shows n(rj, 9b) calculated numerically for point 1 in upwind (curve 1) 
and for point 2 in downwind (curve 2). It is seen that the most part of the 
interstellar atoms arrive to the point 1 from 9 ~ 15° and to the point 2 from 
9 ~ 55°. Therefore, contrary to the simple consideration the most of H atoms 
arrive to the vicinity of the Sun not exactly from upwind. 

Now, it is clear that the excess in the backscattered Ly-a intensity in mod- 
els 1 and 2 as compared with the model 3 is connected with larger number 
densities of H atoms inside the heliosphere in models 1 and 2 due to non vary- 
ing with 9 boundary conditions while in model 3 the number density at 90 
AU decreases with 9. 

It also becomes clear why the intensities (and line-shifts) of the radiation 
scattered by the primary interstellar atoms are in agreement for all models 
(Fig. [2]B). This is because on the one hand, angle-dependence of hydrogen 
parameters at 90 AU for the primary interstellar atoms is weaker than for the 
secondary interstellar atoms, and on the other hand, distribution function of 
primary interstellar atoms is closer to Maxwellian than distribution function 
of secondary interstellar atoms. 

From the middle column of Fig. |2] (A,C) it is seen that there are noticeable 
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differences in the line-shifts of models 2 and 3 in the downwind direction. 
The differences are seen for the total radiation as well as for the radiation 
scattered by the secondary interstellar atoms. At the same time there is almost 
no discrepancies in line-shifts of model 1 and model 3. In spite of model 1 is 
the simplest model without any effects of the heliospheric interface, and one 
could expect the difference. 

In order to understand these results we computed the contribution to the 
total radial velocity of hydrogen at point 2 (7 AU in downwind) from the 
particles that reach this point from the different directions. Similar to n(rj, Ob) 
contribution to the radial velocity can be calculated as follows: 

K-Oi, b ) = +QO 1 r-j— • / f{r h Wi) W i;r dWi 

Fig. H]B shows K( r i, Ob), calculated at point 2 for different models. It is seen 
that the maximum of V r (ri, Ob) is located approximately at #=45° for all mod- 
els. However, model 2 (dashed curve in Fig. 3 B) has a non-negligible contri- 
bution of the negative values of V r {jCi, Ob) for 9 > 100°. The negative values of 
V r (Yi, Ob) are due to relatively hotter secondary interstellar atoms in the down- 
wind hemisphere (due to large temperature of secondary interstellar atoms and 
absence of decrease of the number density from upwind to downwind in model 
2), which can reach point 2 from large values of angle 0. This makes the total 
line-shift of model 2 in the downwind region smaller as compared with models 
1 and 3. Contributions of negative V r (ri, 9b) in model 1 and model 3 are smaller 
than in model 2 but due to different reasons. In model 3 the contribution of 
the particles with negative K( r «, Ob) are significantly reduced because of small 
number density of such particles at the outer boundary that follows from the 
self-consistent model results. As for model 1 we do not see this effect due to 
relatively smaller temperature of the mixture of primary and secondary inter- 
stellar atoms as compared with the temperature of the secondary interstellar 
atoms which exist in model 2. 

Now let us consider the line- width of the Ly-a (right column in Fig. [2]). 
Plot A demonstrates that there are essential qualitative differences in the line- 
widths calculated on the basis of the three models. Model 1 shows a monotonic 
increase of the line-width with the angle 0. Model 2 shows a minimum of the 
line- width at ~60°. Model 3 shows a small local maximum at ~150°. These 
results demonstrate that the line-width of the backscattered Ly-ct profiles is 
very sensitive to the properties of hydrogen distribution at the entrance to 
the heliosphere. In the next section we will focus on understanding of these 
qualitative differences and will show that the main causes of the difference 
are in kinetic non-Maxwellian nature of the hydrogen velocity distribution 
function at 90 AU. Note also, that for the line-widths the simplified model 
2a results are close to model 3 in the upwind hemisphere. However, large 
discrepancies in the downwind region still take place. 
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Now we will consider the differences in the line- widths of the backscattered 
spectra calculated separately for the primary and secondary populations of H 
atoms (right columns in Fig. [2]B,C). For both populations we see the large 
maxima of the line- widths at 9 ~ 90° in the results of model 3. Such strong 
maxima do not exist for model 2 and for model 2a. The latter means that the 
maxima are not due to ^-dependence of hydrogen parameters at 90 AU but 
due to non-Maxwellian veloci ty distribution f unctio n at the outer boundary. 



Namely, as it was shown by iBaranov et al.l (119981 ) (see also, Izmodenov et 



a!., 2001, and Fig. 3 e,f in Katushkina and Izmodenov, 2010) components T 2 
and T p of the kinetic temperatures calculated from the velocity distribution 
function (i.e. T 2 (r) ~ J /(r, w)-(V z (r)-w z ) 2 dw and T p (r) ~ / f(r, w)-(V p (r)- 
w p ) 2 dw) are essentially different. In other words the mean thermal velocities 
of the H atoms are different in different directions. Moreover this difference 
between T 2 and T p temperatures increases approaching to the Sun due to 
local effects. The large maxima of the line-widths at 9 = 90° (Fig. [2]B,C) are 
connected with the changes of the radial kinetic temperature T r of H atoms 
at 90 AU as a function of 9. For example, when observer looks toward upwind 
(i.e. 9 = 0°) T r =T 2 , for the line-of-sight of 9 = 90° T r =T p and for line-of- 
sight of 9 = 180° T, r =T 2 again. In model 3 at the outer boundary T p is larger 
than T z for each of the interstellar populations of H atoms (that follows from 
Baranov-Malama model results). This leads to the maxima of T r at 9 = 90° 
and that are reflected in the maxima of the line-widths for photons scattered 
by each population of interstellar H atoms separately. 



Finally, to qualify the differences in the Ly-a intensities, line-shifts and 
line-widths calculated by making use of models 1-3 we present Table [T] that 
shows differences (absolute and relative) between the model results for the 
upwind (9 = 0°), crosswind (9 = 90°) and downwind (9 = 180°) directions. 
It is seen that one-component hot model 1 gives from 2 % (in upwind) to 
24 % (in downwind) differences with model 3 in intensities and about 10-18 % 
differences in the line-widths. Two-component hot model 2 leads to 8-27 % 
discrepancies with model 3 in the intensities, and from 10 % (in upwind) to 
40 % (in downwind) discrepancies in the line-widths. The differences in the 
line-shifts of models 1 and 2 relative to model 3 are large especially for the 
crosswind direction. However, they are connected with the small values of the 
line-shifts and most probably can not be detected experimentally. 



In the next section we discuss nature of the dependence of the Ly-a line- 
widths on the angle 9. 
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Table 1 

The comparison of intensities, line-shifts and line- widths of the backscattered Ly- 
a radiation for three considered models of the hydrogen distribution inside the 
heliosphere (description of the models see in text) 



6 = 0° 



Ilos{R) 



Vi os (km/s) 



Ties (K) 



lItOS r Il ° S ' 31 • 100% 



\Vlos — Vios, 3l 



100% 



100% 



model 1 
model 2 
model 3 



677 
720 
662 



-24.2 
-24.5 
-24.5 



14 710 
14 158 
12 841 



2.3 
8.7 
0.0 



1.2 
0.0 
0.0 



14.5 
10.3 
0.0 



90° 



Ilos(R) 



Vios (km/s) 



Tios (K) 



I ^los ^los.Z I 



• 100% 



Vl os 3 



• 100% 



• 100% 



model 1 
model 2 
model 3 



398 
421 
380 



-0.9 
-0.6 
-0.3 



17 505 
15 561 
14 834 



4.7 
10.8 
0.0 



200.0 
100.0 
0.0 



18.0 
4.9 
0.0 



180° 



hos{R) 



Vios (km/s) 



Tios (K) 



I ^los ^los,3 I 



100% 



I Vlos — Vlos, 3l 

Vina 3 



100% 



100% 



model 1 
model 2 
model 3 



127 
129 
102 



20.3 
17.4 
20.6 



18 984 
23 946 
17 060 



24.5 
26.5 
0.0 



1.4 
15.5 
0.0 



11.3 
40.4 
0.0 



5 Line-width of the backscattered Ly-cc profiles as diagnostics of 
the heliospheric interface nature 



Costa et al.l ( 119991 ) and lQuemerais et al.l ( 120061 ) have analyzed spectral prop- 
erties of the backscattered solar Ly-o radiation measured by SOHO/SWAN in 
1996-2003. The Ly-o line-width as function of the line-of-sight direction was 
studied. It was shown that there is a noticeable minimum of the line-width 
at 9 — 50° — 60°. This minimum was interpreted as the indication of the two 
(primary and secondary) components of the interstellar H atoms in the helio- 
sphere and, therefore, as an evidence of the heliospheric interface. Indeed, the 
results obtained in the frame of the one-component classical hot model show 
(dashed-dot curve on Fig. 2A) the monotonic increase of the line-width from 
upwind to downwind. Existence of the two components, which have the bulk 
velocities shifted one with respect to other, can help to produce the minimum 
of the temperature in the directions close to the crosswind direction. 

Models 2 and 3 considered in this paper have the two-components of H 
atoms at the outer boundary at 90 AU. Therefore, one can expect that the 
minimum of the line-widths should be obtained in these models. However, as 
it is seen from Fig. [2] A the minimum is seen only for model 2 and not seen 
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i — 1 — r 

120 150 



Figure 5. Results of test calculations of the Ly-a line-width as a function of the 
angle 8. Results based on the one-component hot model (model 1 curve 1) and 
the two-component hot model (model 2 curves 2 and 2b) are shown in plot (A) 
and results based on model 3 (curves 3, 3b, 3c) are shown in plot (B). Curves 2 
and 2b correspond to calculations based on model 2 with and without ionization, 
respectively; curve 3 corresponds to model 3, and curve 3b and 3c are the results of 
model 3 where the i onization was switche d off; in all calculations except model 3c 
solar spectrum from Lemaire et al. ( 19981 ) was assumed, for model 3c the flat solar 
spectrum was employed. Also see Table [2l 



for model 3. Instead, for model 3 we see a small maximum of the temperature 
at 9 = 150°. In order to explain these features we have performed series of 
additional test calculations. 



Our goal is to understand and distinguish roles of different effects influenc- 
ing the line-width of the Ly-a profiles. To do that we studied the effects of all 
possible model parameters and established that the following factors influence 
the dependence of the line-width (T/ os ) on the angle 9: 

1. Ionization processes that change the parameters of H atoms near the Sun. 
From our test calculations we found that the solar gravitation and radiation 
pressure have smaller influence on the line-width than ionization. 

2. Existence of the two populations of the interstellar hydrogen atoms that 
are shifted in velocity space. This effect leads to appearance of the minimum 
of Ti os at 9 = 90° as it was discussed before. 

3. Non-Maxwellian behavior of the velocity distribution function of the two 
populations at 90 AU. More precisely, the large difference between T z and 
T p "temperature" components (z and p are cylindrical coordinates). This 
effect takes place for model 3. 

4. Shape of the solar spectrum of the Ly-a radiation. It will be shown below 
that the shape of the solar spectrum has some small but interesting effect 
on the line-width as function of 9. 
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Fig. [5] summaries the results of test-calculations (see, also, table|2]for model 
numeration) and explore the listed above effects. For all calculations of this 
section we assume the solar gravitation force to be in balance with the solar 
radiation, i.e. the parameter /i = 1 and the trajectories of the H atoms are 
straight lines. To explore the effects of ionization in the vicinity of the Sun we 
performed test calculations with typical ionization rate of (3e 
as well as the calculations with /3g = (models 2b, 3b, 3c). 

Table 2 

Parameters of the test calculations 



5.9 • 10" 7 s 



curve's number 
on fig. [5] 


Model description 


1 


one-component model (1) 


2 


two-component model (2) 


2b 


two-component model (2) without ionization 


3 


our model (3) 


3b 


our model (3) without ionization 


3c 


our model (3) without ionization; flat solar spectrum of the Ly-o radiation 



The effects of the local ionization on the function of Ti os (9) are clearly il- 
lustrated in the frame of the one-component hot model 1 (curve 1 in Fig. [5] A) 
since neither two-component nor non-Maxwellian effects are taken into account 
in this model. The line-width increases monotonically from the upwind direc- 
tion to downwind. Such a behavior of the line- width of the backscattered Ly-a 
profile reflects the properties of the H atom distribution since the ionization 
(due to so-called selection effect) leads to decrease of effective temperature of 
H atoms in upwind and incr ease of the temper ature in the downwind direction 
(see for example Fig. 3.5 in llzmodenovl (120061 )). 



Results of the two-component hot model 2b that does not account for 
ionization (curve 2b on Fig. [5]) illustrate the effect of the two populations 
of the interstellar H atoms. These two populations have different bulk ve- 
locities V g (at 90 AU in the upwind direction: V ZjPrimary ~ —27 km/s and 
Vz,secondary ~ — 16 km/s) and rather small thermal velocities. Therefore in 
the upwind and downwind directions (where V r = ±V Z ) the line-of-sight pro- 
jections of the velocity distribution functions of the primary and secondary 
components are overlap each other in the velocity space only partially. In the 
crosswind direction V r = V p ~ for both primary and secondary interstellar 
atoms. This means that in this direction the projections of the distribution 
functions on the radial line-of-sight overlap completely. That is why the ra- 
dial temperature of the mixture of primary and secondary interstellar atoms 
is smaller in the crosswind direction than in the upwind and downwind di- 
rections. This minimum of the radial temperature of H atoms in crosswind is 
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reflected in the Ly-a line-widths as it is clearly seen from model 2b. 



Results of the two-component hot model (curve 2) where the ionization is 
taken into account combine both the increase of T[ os from upwind to down- 
wind due to ionization effect and the local minimum in crosswind due to two 
populations of H atoms. Hence, a small minimum of TJ 0S at 8=50-60° is seen 
in curve 2. 

The line-width obtained by making use of model 3 is presented in Fig. 0B. 
Note that this model takes into account all considered effects of the helio- 
spheric interface, namely: 1) two populations of interstellar H atoms, 2) 9- 
dependence of the hydrogen parameters at 90 AU, 3) non-Maxwellian features 
of the hydrogen velocity distribution function at the entrance to the helio- 
sphere. In particular, differences between kinetic temperature's components 
T z and T p play an important role. 

Curve 3c corresponds to the model 3c that does not take into account 
ionization and in addition in model 3c we assume that the solar Ly-a flux 
does not depend on freque ncy. Note that i n all o ther models we use the shape 



of the solar spectrum from lLemaire et al.l (119981 ) 



It is seen that the effect of minimum in 7} os at 9 = 90° (due to two popu- 
lations) practically disappears in model 3c, but still is visible. This practical 
disappearance of the minimum is connected with the increase of the radial ki- 
netic temperature T r with increase 9 at 90 AU. The effect of T r increase with 
9 compensates the effect of the minimum of the radial hydrogen temperature 
in crosswind due to the two populations of H atoms. That is why the value 
of the minimum of the line-widths at 9 = 90° is much smaller for model 3c as 
compared with model 2b. 

It is interesting to note, that for model 3b (that is more realistic than model 
3c) local minimum of 7] os at 9 = 90° replaces by small maximum. This effect 
is due to shape of the solar spectrum. Remember, that in model 3c we use 
flat solar spec t rum w hile in model 3b we use nonuniform solar spectra from 



Lemaire et al.l (119981 ) . 



Let us return to model 3 without any additional assumptions. It is seen 
(curve 3) that there is no minimum of Ti os at 50-60° at all in this case. But 
there is a small maximum at 9 = 150°. This maximum can be explained 
by ^-dependence of the kinetic temperature T r . The radial temperature of 
sum of the primary and secondary interstellar atoms is shown in Fig. [H] (plot 
A corresponds to model 2, plot B - to model 3). It is seen that in case of 
the two-component hot model 2 ionization leads to maximum of the T r in 
downwind. For model 3 the maximum of T r is located at about 9 = 150°. This 
effect is reflected in the line-widths of Ly-a radiation that is seen in curve 3 
of Fig.|BjB. 
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z, AU 



z,AU 



Figure 6. The kinetic radial temperature (T r ) for the sum of primary and secondary 
interstellar atoms in the heliosphere; Plot A corresponds to the two-component hot 
model 2, plot B corresponds to the results of the model 3. In these calculations 
p b = \^ E = 5.9 . 10-7S- 1 . 

6 Conclusions and discussion 



In this work we applied different models of hydrogen distribution in the 
heliosphere to compute the spectral properties of the backscattered solar Ly-cv 
radiation as it would be measured at 1 AU in the anti-solar directions. We have 
found out that imprints of the heliospheric interface in the H atom distribution 
inside the heliosphere have a strong influence on the Ly-a parameters. 



We considered the minimum of the line-width of the backsc attered Ly-cy 



radiation at 50-60° from upw ind that was observed by SWAN (j Cost a et al. 



19991 ; iQuemerais et all 120061 ). In the experimental data the line- width in the 
directions of 9 =50-60° is smaller than in the upwind direction by 1500-2000 K. 
This minimum is seen for 1996 and 1997 and practically is not seen for 2002- 



Quemerais et al.. 


2006 


1999 




Quemerais 


et al. 



. Fig. 7 ) . This minimum was explained (see I Cost a et al. 



20061 ) by existence of the two different populations of 



the interstellar hydrogen atoms that are shifted in velocity space. However, we 
noticed that the line- width calculated with 2D stationary Baranov-Malama 
model (IQuemerais fc Izmodenovi . 120021) has no any minimum of T; os . Non- 
stationary 2D Baranov-Malama model (j Quemerais et all 120081 ) predicts small 
minimum in 2003, but there is no minimum in 1997. 



In this work we theoretically explore the nature of the observed minimum 
on the basis of three types of the hydrogen distribution inside the heliosphere. 
It was shown that the minimum of the line-widths appears only for the two- 
component hot model, and there is no minimum at all for model 3 that takes 
into account all effects of the heliospheric interface. It was found that the 
absence of the minimum in model 3 is because the effect of the two components 
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is compensated by the non-Maxwellian features of the velocity distribution of 
H atoms at 90 AU after they passed the heliospheric interface region, namely, 
by strong anisotropy of the kinetic temperatures of H atoms (T z < T p ). 

Therefore, the question why the the minimum of the Ly-a line-width exists 
in the experimental data remains to be open. Possibilities to get the mini- 
mum still exist in the frame of model 3. Firstly, models considered here does 
not take into account the effects of latitudinal and solar cycle variations of 
the photoionzation and charge exchange rates as well as the solar radiation 
pressure. These local effects may potentially change the result of this paper. 



Another possibility is to change the boundary conditions at 90 AU, i.e. to 
make a change in the model of the heliosp heric interface. For exam ple, inter- 
stellar magnetic field may play a key role (jlzmodenov et al.l . 120051 ) . Another 



possibility is the multi-compon e nt na t ure of both the h e liospheric and inter- 



stella r plasmas ( IMalama et al.l . 120061 ; llzmodenov et al.l . l2009t IChalov et al 



20101 ) . In this non-equilibrium plasma model the interstellar pickup ions are 
treated as the separate kinetic component. The plasma temperature in the 
vicinity of the heliopause is smaller in this model as compared with the 
Baranov-Malama model. Therefore, we could expect a decrease in the kinetic 
temperatures of the secondary interstellar atoms. This might make larger ve- 
locity space separation of the primary and secondary interstellar atoms at 
90 AU. The separation may effectively increase the effect of the two popula- 
tions as compared with the other effects. As the result one may hope that the 
observed minimum will appear in the model. 



All these possibilities will be checked in the nearest future. 

This work was partially supported by RFBR grants 10-02-93113, 10-02- 
01316 and "Dynastia" foundation. Part of this work was done in the frame of 
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